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Abstract:  The  study  aims  are  the  active  control  of 
transonic  buffet  over  airfoils  and  wings.  A new  mov- 
ing device  at  the  trailing-edge,  a so-called  “Trailing 
Edge  Deflector”  (TED)  designed  at  ONERA,  is  used 
to  decrease  the  buffet  instabilities. 

This  study  is  limited  to  buffet  phenomenon  on  upper 
surface  of  a stiff  2D  airfoil  in  transonic  flow. 
Experimental  and  numerical  results  allow  to  better 
understand  the  phenomena  and  to  validate  the 
efficiency  of  the  TED. 

The  TED  can  be  moved  sequentially  to  a static  posi- 
tion, or  better  be  driven  by  dynamic  servo- motions. 
A selected  static  position  increases  aerodynamic 
performances  for  high  lift  coefficient  and  delays  the 
buffet  onset  vs.  lift.  A preferable  TED  active  control 
law,  coupled  with  measured  flow  signals,  allows  to 
greatly  reduce  the  aerodynamic  instabilities. 

The  time-consistent  strong  viscous-inviscia  in- 
teraction approach  VIS  15,  developed  previously 
for  computing  buffet  or  unsteady  separated  flows 
over  airfoils,  gives  access  to  a direct  time-accurate 
simulation  of  the  active  buffet  control  by  TED 
servo-command,  both  in  open-loop  control,  and 
in  the  closed-loop  servo-command  driven  by  the 
computed  self-induced  signals  of  buffet. 

Key  words:  Unsteady  flow,  buffet,  buffet  control, 

two-dimensional  flow , shock-wave  interaction,  flow  sep- 
aration, numerical  method,  strong  interaction , viscous- 
inviscid  interaction,  active  control,  control  law,  closed- 
loop  control,  open-loop  control . 

List  of  Symbols 

c airfoil  chord,  c = 200  mm 

x longitudinal  coordinate 

2 vertical  coordinate 

M Mach  number 

U free-stream  velocity 

T temperature 

P pressure 

Cr,  pressure  coefficient 

Kec,c  Reynolds  number 

p density 

t time 

/,  f0  frequency,  natural  buffet  frequency 
<j>  phase 

Xsh  shock  wave  position 

a angle  of  attack  (degrees) 

8 angle  of  deflector  (degrees) 

8mean  static  position  of  the  deflector 

A amplification  of  the  control  law 

t time  delay  of  the  control  law 


Subscripts 

t.)0  free  stream  value 

. h total  value 

. )def  deflector  value 

J mean  mean  value 

| . [ amplitude  function  for  fluctuations 

1.  Introduction 

Buffeting  is  characterized  by  structure  vibrations 
induced  by  aerodynamic  excitations.  The  buffeting 
phenomenon  can  exist  on  different  parts  of  the  flows 
around  vehicles.  It  can  appear  on  aircraft,  rockets, 
turbomachine  stages...  In  this  paper,  the  study  is 
focused  on  the  upper  surface  of  a wing  in  a transonic 
flow. 

The  shock  wave/ turbulent  boundary  layer  interac- 
tion and  flow  separations  induce  flow  instabilities, 
named  buffet.  It  results  in  pressure  fluctuations 
on  the  wing,  that  create  unsteady  loads  which 
can  cause  structure  vibrations  on  its  eigen  modes. 
The  buffet  loads  may  cause  fatigue  problems.  In 
addition,  maneuverability  and  handling  qualities  are 
often  degraded. 

The  buffet  phenomenon  appears  when  the  Mach 
number  or  the  angle  of  attack  of  the  aircraft 
increases.  The  flight  envelope  is  limited  by  this 
phenomenon. 

Of  course,  any  delay  of  the  buffet  onset  by  de- 
creasing or  canceling  the  aerodynamic  instabilities 
would  permit  to  increase  flight  performances  of  the 
aircraft.  For  example,  the  aircraft  take-off  weight, 
maneuverability,  or  the  range,  can  be  increased. 

The  aim  of  this  study  is  to  control  the  buffet 
phenomenon  by  a direct  action  device  on  the  flow 
instability.  As  a new  moving  part  at  the  trailing  edge 
of  the  wing,  a so-called  “Trailing  Edge  Deflector” 
(TED),  designed  at  ONERA,  is  used  to  decrease 
these  instabilities. 

This  control  system  study,  here  limited  to  two- 
dimensional  flows  and  to  the  aerodynamic  aspect 
of  the  phenomenon,  will  be  further  studied  in 
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three-dimensional  flow  [1,  2]. 

In  order  to  analyze  the  effect  of  the  new  control 
system  on  the  buffet  instabilities  and  to  help  for 
better  understanding  of  the  manner  in  which  the 
deflector  can  act  on  flow  separations  and  shock 
wave  oscillations  over  the  airfoil,  experimental 
investigations  were  performed  on  the  stiff  2D  airfoil 
(OAT15A)  in  the  T2  wind  tunnel. 

Numerically  in  addition,  the  time-consistent  strong 
viscous-inviscid  interaction  approach  VIS  15  [3,  4], 
developed  previously  for  computing  buffet  or  un- 
steady separated  flows  over  airfoils,  was  used  to 
give  access  to  a direct  time-accurate  simulation  of 
the  active  buffet  control  with  TED  servo-command, 
both  in  open-loop  control,  and  in  the  closed-loop 
servo-command  driven  by  the  computed  self-induced 
signals  of  buffet  themselves  [5]. 

2.  Model  and  instrumentation 

An  experimental  investigation  on  the  effects  of  the 
deflector  on  buffet  characteristics  of  an  OAT  15 A 
airfoil  was  carried  out  in  T2  wind  tunnel  of  DMAE, 
a department  at  ONERA  Toulouse  (fig:l).  T2  is 
a transonic,  pressurized  and  cryogenic  wind  tunnel 
with  closed  circuit  [6].  The  test  section,  equipped 
with  top  and  bottom  flexible  walls  is  0,37m  high  by 
0,39m  wide  and  1,4m  long.  The  OAT15A  airfoil 
is  a supercritical  airfoil  designed  by  ONERA  with 
a thickness-to-chord  ratio  of  12,6%  , a chord  of 
200mm  and  a span  of  0,39m  extending  the  full 
width  of  the  tunnel.  The  model  has  a design  point 
of  M — 0, 73  and  a lift  coefficient  of  Ci  = 0, 65.  It  is 
equipped  with  60  steady  and  19  unsteady  pressure 
transducers.  Their  locations  on  the  airfoil  are  shown 
in  figure  2.  The  signals  were  sampled  at  15kHz  with 
a filter  at  bkHz.  All  of  the  tests  were  performed  at 
a Reynolds  number  of  approximately  Rec  = 4,  5.106 
with  transition  fixed  at  | = 7%  on  upper  and  lower 
sides  of  the  model.  The  tests  were  performed  at 
ambient  temperature.  The  Mach  number  varied 
between  M0  = 0, 72  — 0, 78  and  the  maximum  value 
of  the  angle  of  incidence  was  a = 3,5°.  For  the 
design  Mach  number  of  M0  = 0,  73,  the  shock  wave 
position  moves  with  a frequency  of  around  / = 75 
Hz  in  “buffet”  conditions. 

• Trailing-edge  deflector  description 

The  deflector  is  situated  at  the  lower  surface  of  the 
trailing  edge  (fig:3).  It.  is  maintained  along  the 
wingspan.  Its  chord- wise  length  is  only  1 to  3 percent 
of  the  airfoil  chord.  It  can  be  moved  from  -5°  up  to 
45°,  and  have  a maximum  frequency  of  200Hz.  For 
manufacturing  problems,  it  is  made  of  3 parts.  Each 
of  them,  simultaneously,  are  set  to  the  same  position 
or  put  in  motion.  It  is  driven  by  an  electrical  mo- 


tor. Aerodynamic  results  of  the  deflector  effects  are 
presented  in  the  next  paragraphs. 

3.  Upper  surface  buffet  description  in 
two-dimensional  transonic  flow 

Buffeting  can  appear  in  many  flow  conditions.  In 
transonic  flow,  it  is  accentuated  by  the  shock  wave 
motion. 

Natural  instabilities  on  airfoil  appear  by  an  increase 
in  the  Mach  number  or  the  angle  of  attack.  In  tran- 
sonic flow,  the  buffet  is  characterized  by  a complex 
interaction  between  shock  wave  and  boundary  layer 
[7,  8].  The  shock  wave  oscillations  are  caused  by 
flow  separation,  as  shown  by  the  numerical  solution 
VIS  15  (fig:4).  The  buffet  onset  starts  when  the  sep- 
aration bubble  increases  and  joins  the  trailing  edge. 
Then,  instabilities  develop  on  a large  scale. 

The  shock-wave  moves  upstream,  pushed  by  flow  sep- 
aration. Therefore,  the  flow  separation  decreases  and 
the  shock  wave  moves  back  downstream.  Thus,  the 
intensity  of  the  shock  wave  increases  and  a new  sepa- 
ration bubble  is  created.  The  flow  separation  spreads 
to  the  trailing  edge  and  push  the  shock  wave  up- 
stream again  (fig: 4). 

One  of  the  characteristics  of  the  buffet  in  two- 
dimensional  flow  is  that  the  harmonic  spectrum  of 
the  shock  wave  motion  exhibits  a dominant  fre- 
quency, plus  lower  harmonics.  This  fundamental  fre- 
quency depends  on  the  shape  of  the  airfoil  and  on 
the  aerodynamic  flow  conditions  [9,  10]. 

4.  Time-accurate  numerical  simulation 
of  active  buffet  control  : VIS  15  code 

The  VIS  15  numerical  method  is  a “time-consistent 
strong”  Viscous-inviscid  Interaction  (VII)  approach 
that  was  generated  a few  years  ago  [3,  4],  and  that, 
was  the  first  VII  method  giving  the  same  capability 
of  direct  access  to  a time-accurate  simulation  of  tran- 
sonic buffet  [3,  11,  12]  or  unsteady  separated  flows 
over  airfoils  [3,  4]  as  the  RANS  approaches.  More 
recently,  the  VIS  15  method  was  extended  to  the  EU- 
ROSHOCK passive/active  control  of  buffet  through 
porous  walls  [13],  and  to  the  full  time-accurate  direct 
simulation  of  active  buffet  control  in  closed  servo- 
command  loop  with  the  present  TED  device  [5],  or 
with  the  ERE  A deforming  bump  device. 

The  VIS  15  code  is  basically  issued  from  the  common 
methodology  of  the  Viscous-inviscid  Solvers  “VIS” 
developed  at  ONERA  from  1979  [14],  using  in  un- 
steady equations  the  same  turbulent  models  and  vis- 
cous numerical  schemes  as  first  designed  in  steady 
2D-3D  viscous  flows  [11,  15,  16,  17],  but  restrict- 
ing (in  VIS  15)  the  in  viscid  part  of  the  VIS  solver 
to  the  transonic  small  perturbation  approximation, 
while  extending  the  VII  strong  coupling  to  a time- 
consistent  one  (viscous  and  inviscid  unsteady  equa- 
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tions  “discretized  and  coupled  until  convergence”  at 
same  time-step),  so  recovering  the  same  domain  of 
influence  and  same  fully  time-parabolic  nature  as  the 
averaged  Navier-Stokes  RANS  equations. 

In  VIS  15,  a thin-layer  approximation  of  Le  Balleur’s 
“Defect-Formulation  theory”  [14,  15,  17,  16]  for 
Navier-Stokes  equations  is  used.  This  theory  as- 
sumes at  all  points  a splitting  of  the  equations  sys- 
tem (RANS  Navier-Stokes)  in  two  exactly  equiva- 
lent strongly  interacting  equation  systems  : the  vis- 
cous “Defect-Formulation”  system  and  the  pseudo- 
inviscid-flow,  both  of  them  being  solved  on  the  same 
physical  domain.  The  advantages  of  this  VII  ap- 
proach consist  in  a very  low  computational  cost  that 
doesn’t  grow  with  the  Reynolds’  number,  a much 
lower  undesiderable  numerical  viscosity,  and  a high 
flexibility  in  choosing  at  will  the  viscous  approxi- 
mations level,  ranging  from  boundary-layer  to  full 
RANS.  At  any  approximation  level,  the  full  overlay 
of  viscous  and  inviscid  computational  domains  en- 
sures the  respect  of  the  full  viscous  upstream  influ- 
ence in  supersonic  areas  [18],  and  of  a normal  pres- 
sure gradient  inside  viscous  layers  [14,  15,  17,  16]. 
The  thin-layer  viscous  approximation  in  VIS  1 5 is 
assumed  in  a “displacement-body  reference  frame” 
[17,  16],  which  extends  the  range  of  the  thin-layer  as- 
sumption and  gives  access  to  “massive  separations” 
and  “deep  stall” . 

The  viscous  numerical  method  [14,  16,  17,  3]  is  a 
field/integral  hybrid  method,  with  y-discretization 
(but  with  y-rank  lower  than  discretization),  that  in- 
volves a space-marching  technique  using  non-linearly 
implicit  schemes,  implicit  turbulence  coupling,  and  a 
switch  of  direct/inverse  modes  of  solution.  It  is  us- 
able both  in  steady  or  unsteady  conditions.  At  each 
viscous  station,  “parametric”  modeled  turbulent  ve- 
locity profiles  [17,  16,  14,  15]  — ranging  from  at- 
tached flow  to  massive  separation  — are  discretized 
along  the  direction  normal  to  the  local  inviscid  in- 
teracting flow  streamlines,  on  a normal  grid  that  is 
self-adaptive  to  the  boundary  layer  thickness  and  the 
maximal  normal  velocity  gradient  (fig.9). 

The  method  includes  first  an  algebraic  turbulence 
model  [14,  15,  17,  16],  based  on  the  discretized  para- 
metric velocity  profiles  and  on  a mixing- length,  and 
secondly  an  out-of-equilibrium  two  equation  model 
“k  - uV  forced  ”,  Le  Balleur  [14,  15,  17,  16], 
“forced”  by  the  parametric  velocity  profiles  mod- 
elling, and  selected  in  all  present  computations. 

For  the  resolution  of  the  inviscid  system,  an  unsteady 
transonic  small  perturbations  (TSP)  method,  exten- 
sion of  [19]  with  revised  schemes,  is  used  in  VIS15. 
The  strong  coupling  is  numerically  obtained  using 
the  original  “Semi- Implicit”  algorithm  defined  in 
1984  [3,  4],  which  converges  the  VII  coupling  resid- 
uals both  in  direct  and  inverse  viscous  modes,  after 
extracting  a viscous  influence  operator,  and  with  an 


iterative  marching  conditioning.  The  convergence  is 
imposed  at  each  time  step  in  order  to  obtain  a time- 
consistent  strong-coupling  method  and  thus  to  take 
into  account  at  each  time  step  the  full  downstream- 
upstream  viscous  influence  (as  in  Navier-Stokes  so- 
lutions). Doing  this  enables  the  proper  treatment 
of  strong  viscous  interaction  phenomena  in  the  un- 
steady case,  such  as  unsteady  separations  and  shock 
- boundary  layer  interactions. 

5.  Comparison  of  numerical  solutions 
with  experiments 

5.1.  Numerical  conditions 

Numerical  flow  conditions  are  identical  to  experi- 
ments (table  1). 


Rcc 

Pto 

Tto 

c 

xtr/c  | 

4,3  106 

1.7  bar 

300  K 

0.2  m 

7% 

Table  1.  Numerical  conditions 


The  VIS  15  inviscid  grid-part  used  for  calculations 
is  (200  x 100),  with  a few  streamwise  refinements 
at  leading-edge,  trailing-edge,  and  in  the  region 
of  oscillating  shock-induced  separation.  Figure  8 
shows  a zoom  of  the  inviscid  grid.  The  extent 
of  the  computational  domain  is  10  chord  lengths 
upstream/downstream,  and  10  chord  lengths  up- 
per/lower. The  VIS  15  overlaying  viscous  grid-part, 
figure  9,  is  made  of  2 x ( 163  x 49  ) points,  has 
wall  nodes  in  coincidence  with  the  2 x (163) 
inviscid  nodes  along  the  upper /lower  sides  of  the 
airfoil  and  wake-cut,  and  seif-adapts  its  49  nodes  in 
normal  direction  to  the  local  thickness  and  maximal 
gradient  of  the  viscous  velocity  profiles.  The  VIS  15 
grid  refinement  at  shock  separation  is  designed, 
with  local  streamwise  meshes  smaller  than  the 
boundary  layer  thickness  (AX  about  0.5<5),  in  order 
to  effectively  discretize  and  resolve  numerically  the 
instantaneous  physical  extent  of  compressions  at 
turbulent  shock-induced  separation. 

To  allow  a good  comparison  between  experiments 
and  calculations,  it  is  important  to  apply  some 
corrections  on  the  Mach  number  M0  and  on  the 
angle  of  attack. 

T2  wrind  tunnel  is  equipped  wflth  top  and  bottom 
flexible  walls.  As  a result,  classical  corrections  for 
this  two  walls  aren’t  necessary,  at  least  at  steady 
flow  conditions.  However  boundary  layers  that 
appear  on  the  lateral  walls  of  the  wind  tunnel  induce 
a flow  deceleration.  A negative  correction  of  the 
Mach  number  is  necessary.  For  T2  wind  tunnel,  this 
correction  is  estimated  to  5M  = -0.06. 

For  a good  comparison  it  is  essential  to  compare  the 
same  level  of  buffet.  So  a correction  of  the  angle 
of  attack  is  done  for  the  two  particular  studied 
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cases  : one  before  the  buffet  onset  and  one  another 
beyond  buffet  onset.  Corrections  applied  to  this  two 
cases,  which  are  still  not  fully  optimized  (in  shock 
position),  are  indicated  in  table  2. 


case 

iff  Oexpe 

MoviS  15 

OLexpe. 

OV/S’15 

1 

0.736 

0.73 

2.5° 

2.34° 

2 

0.736 

0.73 

2.125° 

2.1° 

Table  2.  Studied  cases 


5.2.  Characteristics  of  OAT15A  buffet 

The  transonic  buffet  is  characterized  by  pressure 
fluctuations  and  oscillations  of  the  shock  wave. 
The  figure  5 shows  the  pressure  fluctuations  on 
the  airfoil.  The  unsteady  pressure  transducers  on 
the  oscillation  band  of  the  shock  wave  indicate  a 
high  level  of  fluctuations.  Downstream  of  the  shock 
wave,  fluctuations  are  still  important  due  to  flowr 
separation  instabilities.  Upstream  of  the  shock 
wave,  in  the  supersonic  region,  the  fluctuations  are 
low. 

The  numerical  VIS15  solution  reproduces  well  the 
fluctuations  behaviour  (fig:5),  despite  an  optimal 
adjustment  of  the  shock  mean  position  is  not  here 
achieved,  as  shown  by  the  zero-harmonic  of  the  pres- 
sure (fig:6),  wdiich  also  displays  the  good  prediction 
of  the  steep  pressure  evolution  at  trailing-edge  on 
lower-side. 

The  instabilities  are  mainly-harmonic.  Figure  7 
illustrates  the  shock  position  vs  time.  The  shock 
wave  oscillation  frequency  is  about  / = IbHz 
in  experiment,  whereas  a somewhat  lower  fre- 
quency (increasing  with  angle  of  attack  from  about 
/ — 53 Hz  for  incipient  buffet  to  / = 60Hz  for  well 
established  buffet)  is  computed  by  VIS  15.  It  is 
not  yet  fully  clear  whether  this  numerical  shift  in 
fundamental  frequency  is  due  to  the  TSP  inviscid 
truncation  still  used  in  the  case  of  VIS  15,  or  to  the 
form  of  the  TSP  variant  used  [19],  or  to  the  free- air 
solution. 

More  informations  on  buffet  measurements  are  given 
in  references  [20,  21,  22]. 

5.3.  Effects  of  a static  TED  change 

The  trailing  edge  deflector  can  be  moved  to  a 
particular  position.  Before  buffet  onset,  increasing 
static  TED  defection  acts  on  shock  wave  position 
and  on  flow  separation  levels.  The  supersonic  level 
is  increased  and  the  shock  wave  moves  downstream. 
Also,  the  trailing  edge  thickness  is  increased  and 
the  flow  is  modified  in  this  region  (fig:  10).  All 
these  effects  are  well  detected  by  VIS  15  numerical 
calculations  (fig:  11). 


The  buffet  onset  can  be  detected  by  the  pressure 
fluctuation  levels.  The  lift  coefficient  can  be  de- 
termined by  pressure  integration.  The  figure  12 
shows  the  buffet  onset.  The  angle  of  attack  of  the 
buffet  onset  is  well  detected  by  the  VIS  15  solution, 
where  however  the  buffet  onset  is  somewhat  more 
brutal  than  in  the  experimental  results.  The  figures 
12  and  13  indicate  the  effects  of  two  positions  of 
the  deflector  ( 8 = 0°  and  6 = 15°)  on  the  buffet 
onset.  The  incidence  of  buffet  onset  is  decreased 
with  a deflector  angle  increase,  but  the  lift  at  buffet 
onset  is  increased.  Both  effects  are  well  predicted 
numerically  by  VIS  15. 

The  calculations  overestimate  slightly  the  lift  co- 
efficient because  of  their  downstream  shock- wave 
position.  The  figure  13  indicates  the  steady  lift 
coefficient  vs  the  angle  of  attack  for  the  two  po- 
sitions of  the  deflector  (d  = 0°  and  S = 15°).  At 
iso-incidence,  an  increased  deflector  angle  increases 
the  lift  coefficient  for  the  airfoil. 

The  figure  14  represents  the  buffet  onset  by  means 
of  the  fluctuation  amplitude  of  a local  pressure 
coefficient  (x/c  = 0.9).  Experimentally,  the  pressure 
fluctuations  are  calculated  for  two  frequency  bands 
[5  - 4000Hz]  and  [50  - lOOfte].  This  allows  to  know 
the  difference  between  the  separation  flow  noise  and 
the  buffet  phenomenon.  VIS  15  calculations  give 
a buffet  onset  that  is  very  close  to  experimental 
buffet  onset  in  the  [50  — lOOiT^]  frequency  band. 
As  a result,  it  seems  that  the  buffet  phenomenon 
is  well  modelized  by  VIS  15.  An  important  effect 
of  the  deflector  is  to  delay  the  buffet  onset  at  iso 
Mach-Cz.  The  figure  14  illustrates  this  benefit 
displayed  by  experiments,  and  well  reproduced  by 
VIS  15  calculations. 

The  drag  polar  clearly  shows  the  increase  in  aerody- 
namic performances  for  high  lift  coefficient  (fig:  15), 
but,  for  lower  lift  coefficient,  there  is  no  profit. 
VIS  15  gives  the  same  results  (fig:  15). 

5.4.  Open  loop  control  - Effects  of  a dynamic 
TED  motion 

The  deflector  can  be  driven  by  dynamic  movements 
around  a static  position.  Different  deflector  motions 
were  tested  to  try  to  understand  the  effects  of  the 
deflector  motion  on  natural  buffet  (fig:  16). 

• Step  motion  of  TED 

Deflector  step  motions  were  tested.  This  kind  of 
motion  is  very  interesting  because  it  allows  to  study 
and  understand  the  crossing  mechanisms  from  a flow 
state  to  another  one.  Different  cases  were  studied. 
The  figure  17  shows  the  buffet  onset  by  variation 
of  the  deflector  angle.  This  figure  indicates  the 
extremal  upstream  and  downstream  positions  of  the 
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shock-wave.  The  shock-wave  oscillation  starts  more 
brutally  by  VIS  15  calculations  than  by  experiments. 
However,  the  comparison  is  very  good. 

An  interesting  case  is  the  switch  from  a steady 
state  to  a buffet  state.  The  figure  18  indicates  the 
shock- wave  oscillations  for  a step  from  6 = to 
8 = 25°.  VIS15  can  visualize  the  whole  flow  field. 
The  figure  19  shows  the  propagation  in  time  of  the 
Mach  number  perturbation  which  is  created  by  the 
step  at  the  time  t = 0 ms.  The  perturbation  takes 
about  1 ms  to  go  upstream  in  the  subsonic  region 
of  the  upper  surface  and  to  reach  the  shock  wave. 
It  seems  to  be  an  acoustic  propagation.  On  the 
lower  surface  of  the  airfoil,  the  perturbation  takes 
about  2 ms  to  go  to  the  upper  surface.  It’s  also  an 
acoustic  propagation.  The  perturbation  goes  round 
the  leading  edge  and  in  the  supersonic  region.  The 
required  time  to  reach  the  shock  wave  is  about  lms. 
It  is  not  a simple  acoustic  propagation.  Presently, 
we  don’t  know  how  to  explain  it. 

* Sine-shaped  TED  motions  at  steady  flow 
conditions 

Firstly,  sine-shaped  signals  prescribed  to  the  de- 
flector motion  were  tested  at  steady  flow  conditions. 
Figure  20  shows  the  shock  wave  amplitude  oscil- 
lations for  different  frequencies  of  the  deflector 
(experimental  results).  The  deflector  motion  creates 
shock  wave  oscillations.  The  amplitude  of  the  shock 
wave  oscillations  and  its  frequency  depend  on  the 
deflector  motion.  The  deflector  motion  creates 
shock  wave  oscillations  at  the  same  frequency  than 
the  deflector.  This  hanging  up  of  frequency  is 
observed  for  all  frequencies  of  the  deflector.  The 
higher  the  amplitude  of  the  deflector  is,  the  greater 
the  amplitude  of  the  shock  wave  oscillations  will  be. 
The  amplitude  of  the  shock  wave  oscillations  varies 
with  the  frequency  of  the  deflector.  The  highest 
amplitude  of  the  shock  wave  oscillation  is  observed 
for  a frequency  of  the  deflector  which  is  near  the 
natural  buffet  frequency.  All  this  remarks  are  well 
predicted  by  numerical  VIS  15  calculations  in  figure 
21.  Except  for  an  amplitude  deflector  of  A — ±10°, 
the  greater  amplitude  of  the  shock  wave  oscillation 
is  observed  for  a frequency  of  the  deflector  which  is 
lower  than  the  computed  natural  buffet  frequency. 

• Sine-shaped  TED  motion  at  natural  buffet 
flow  conditions 

At  flow  conditions  with  natural  buffet,  the  mo- 
tion of  the  deflector  can,  for  certain  amplitudes  and 
frequencies,  impose  its  frequency  to  shock  wave  and 
flow  separations  oscillations.  The  higher  the  ampli- 
tude of  the  deflector  is,  the  greater  the  influence  field 
will  be.  The  figure  22  illustrates  the  domain  of  in- 
fluence in  frequency  for  3 amplitudes  ( A = ±2.5°; 


A = ±5°;  A = ±10°)  by  calculations  and  experi- 
ments. 

ATS  15  calculation  solutions  give  a slightly  greater  do- 
main of  influence  than  experiments,  but  behaviours 
are  identical. 

For  deflector  motion  frequencies  far  from  the  natural 
buffet  frequency,  the  shock  wave  keeps  its  natural 
frequency,  but  the  amplitude  of  oscillations  are  mod- 
ulated by  the  deflector  motion. 

Experimentally,  the  amplitude  of  the  deflector  A = 
±2.5°  have  just  few  effects  on  the  shock  wave  oscil- 
lations while  calculation  solutions  give  more  effects. 
In  experiments,  the  real  position  of  the  deflector  was 
measured  by  a camera  [20].  For  little  amplitudes  or- 
dered, the  real  deflector  amplitude  is  lower  than  the 
commanded  amplitude.  For  A = ±2.5°,  we  find  a 
real  amplitude  about  A = ±1.0°  [21].  For  the  am- 
plitudes ordered  of  A = ±5°  and  A = ±10°,  the  real 
deflector  amplitude  are  near  the  commanded  ampli- 
tudes. Experimental  and  numerical  behaviours  are 
comparable.  For  these  amplitudes,  the  deflector  mo- 
tion causes  an  increase  in  shock  wave  oscillations. 

5.5.  Closed  loop  command  - Active  buffet 
control 

Some  phase  shifts  of  sine-shaped  open  loop  motions 
of  the  TED  (deflector  angle  phase  vs.  shock  position 
phase)  have  permitted  a brief  stabilization  of  the  un- 
steady buffeting  flow.  But,  it’s  impossible  to  control 
the  buffet  phenomenon  only  by  this  kind  of  motion. 
If  an  open  loop  command  allows  to  understand  the 
effects  of  a deflector  motion,  a closed  loop  of  servo- 
command  is  required  to  stabilize  the  shock  wave  os- 
cillations. It  is  based  on  the  unsteady  measurements 
of  pressure  fluctuations  (fig:25).  We  saw  in  the  pre- 
vious sections  that  a decrease  in  the  deflector  angle 
involves  an  upstream  motion  of  the  shock-wave,  and 
an  increase  in  the  deflector  angle  involves  a down- 
stream shock  wave  motion.  So,  the  idea  is  to  move 
the  deflector  in  order  that  it  is  acting  against  the 
natural  shock  wave  movement. 

The  control  law  for  a measured  signal  P is  described 
by: 

8(t)  — 8mean  4-  J (£),  (1) 

with, 

8\t)=A*(P(t^r)-Pmean),  (2) 

Different  control  laws  were  tested  based  on  the  same 
principle.  Different  measured  signals  can  be  used 
for  control  : pressure  fluctuations  taken  on  a trans- 
ducer which  is  located  on  the  shock  wave  oscillations 
region,  or  shock- wave  position,  or  lift.  The  same  sig- 
nals were  also  enabled  to  drive  the  control  law  when 
this  one  was  introduced  in  ATS  15  [5].  The  control 
law’s  parameters  (A,  r)  weren’t  determined  a pri- 
ori. Different  gains  A and  time  delays  r were  looked 
for.  We  can  find  an  optimal  couple  that  depends  on 
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the  conditions  and  the  chosen  driving  signal.  Dif- 
ferent couples  [A  and  r)  were  tested  both  by  the 
experiments,  and  by  the  numerical  VIS15  calcula- 
tions in  closed-loop  self-command.  Experiments  as 
well  VIS15  calculations  showed  [5]  that  a pressure 
signal  coming  from  the  rear  part  of  the  shock  wave 
oscillating  region  is  the  most  efficient  choice  with  the 
present  control  law,  due  to  it’s  persistent  amplitude 
(i.e.  the  pressure  jump  through  shock  wave). 

By  using  for  the  command  signal  the  unsteady  mea- 
surements of  pressure  fluctuations,  the  buffet  and  the 
shock  wave  oscillations  were  greatly  reduced  in  ex- 
periment (fig.  26  and  27).  By  using  the  same  con- 
trol principle,  the  buffet  can  be  completely  reduced 
(fig:28)  in  the  numerical  VIS  15  simulation.  Calcu- 
lations are  less  expensive  than  experiments.  Then, 
different  new  laws  will  be  tested  by  calculations  in 
order  to  find  new  control  laws  for  three-dimensional 
flow. 

6.  Conclusion 

The  experimental  study  has  shown  the  efficiency  of 
a trailing  edge  deflector  on  two-dimensional  flow. 
The  TED  can  be  moved  to  a particular  position  or 
be  driven  by  dynamic  motions.  A selected  static  po- 
sition increases  aerodynamic  performances  for  high 
lift  coefficient  and  delays  the  buffet  onset.  An  active 
control  law  driven  by  unsteady  measurements  allows 
to  greatly  reduce  the  aerodynamic  instabilities. 
Numerical  solutions  obtained  by  the  time-consistent 
viscous-in  viscid  strong  interaction  code  VIS  15  are 
in  a whole  accordance  with  experiments.  The 
self-command  within  VIS  15  computation  of  the  de- 
flector, driven  in  closed- loop  by  the  computed  buffet 
pressure  signals  themselves,  provides  very  good 
results  for  time-accurate  direct  numerical  simulation 
of  the  active  buffet  control.  Numerical  results  allow 
to  help  to  better  understand  the  phenomena  and  to 
give  a second  validation  of  efficiency  of  the  TED. 
Different  new  laws  will  be  tested  by  calculations  in 
order  to  find  new  control  laws  which  could  be  used 
in  three-dimensional  buffet. 

This  new  control  system  will  soon  be  studied  in 
transonic  three-dimensional  flow.  A similar  model 
extended  to  a transport  aircraft  was  designed  and 
manufactured  with  TED.  Experimental  tests  are 
planned  in  2000. 
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Figure  1.  T2  transonic  wind  tunnel  with  OAT1QA 
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Figure  7.  Shock  wave  position  on  the  airfoil  vs  time  - 
Comparisons  of  numerical  VIS15  solution  (M=0.73;  a = 
2.34°;  <5  = 15°)  with  experiments  (M=0.736;  a = 2.5°; 
S = 15°) 
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Figure  10.  Effects  of  the  static  deflector  position  - Ex- 
perimental Results  (M—0.736;  a = 2.0;  6 = 0°;  5 = 15°) 
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Figure  8.  Inviscid  part  of  VIS  15  mesh  - Zoom  on  partial 
view 


Figure  9.  Viscous  part  of  VIS15  mesh  - Zoom  on  partial 
view 


Figure  11.  Effects  of  the  static  deflector  position  - 
Numerical  VIS15  solutions  (M— 0.73;  a = 2,0;  $ = 0°; 
6 = 15°) 


Figure  12.  Variation  of  the  fluctuations  of  lift  coeffi- 
cient with  the  angle  of  attack  - Comparisons  of  numerical 
VIS  15  solutions  with  experiments 
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Figure  13.  Variation  of  lift  coefficient  with  the  angle  of 
attack  - Comparisons  of  numerical  VIS15  solutions  with 
experiments 


Figure  14.  Buffet  onset  - Pressure  fluctuations  x/c=0.9 
- Comparisons  of  numerical  VIS15  solutions  with  exper- 
iments 
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Figure  15.  Drag  polar  - Comparisons  of  numerical 
VIS  15  solutions  with  experiments 
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Figure  16.  Different  motions  of  the  deflector 


Figure  17.  Variation  of  the  extremal  shock  positions 
with  the  angle  of  the  deflector  - Comparisons  of  numerical 
VIS  15  solutions  (M=0.73;  a = 2.1)  with  experiments 
(M=0.736;  a = 2.13) 


Figure  18.  Variation  of  the  shock  position  vs  time  after 
a step  of  deflector  angle  from  S — 5°  to  £ = 25°  - Com- 
parisons of  numerical  VIS15  solution  (M=0.73;  a = 2.1) 
with  experiments  (M =0.736;  a — 2.13) 
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Figure  19.  Propagation  of  the  perturbation  of  Mach 
number  generated  after  a step  of  deflector  angle  from 
6 — 5°  to  5 = 25°  at  a time  t = 0 ms  - Numerical  VIS15 


Figure  22.  Frequency  domain  of  influence  on  buffet  of 
the  deflector  motion  vs  deflector  amplitude  - Compar- 
isons of  numerical  VIS15  solutions  (M=0.73;  a = 2.34; 
dmean  = 15°)  with  experiments  (M=0.736;  a = 2.5; 

dmean  — 15  ) 


solution  (M=0.73;  a = 2.1) 


Figure  20.  Effects  of  sine-shaped  deflector  motion  be- 
fore the  buffet  onset  - Shock  fluctuations  vs  deflector  fre- 
quency at  3 amplitudes  - Experimental  results  (M=0.736; 
a = 2.13;  dmean  - 15°) 


Figure  23.  Effects  of  sine-shaped  deflector  motion  at 
natural  buffet  conditions  - Shock  fluctuations  vs  fre- 
quency at  3 amplitudes  - Experimental  results  (M=0.736; 
Q:  = 2.5;  dmean  — 15°) 
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Figure  21.  Effects  of  sine-shaped  deflector  motion  be- 
fore the  buffet  onset  - Skock  fluctuations  vs  deflector 
frequency  at  3 amplitudes  - Numerical  VIS  15  results 
(M=0.73;  a = 2.1;  dmean  = 15°) 


Figure  24.  Effects  of  sine-shaped  deflector  motion 
at  natural  buffet  conditions  - Shock  fluctuations  vs 
frequency  at  3 amplitudes  - Numerical  VIS  15  results 
(M=0.73;  a = 2.34;  dmean  = 15°) 
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Figure  25.  Principle  of  the  active  control  law  - Closed- 
loop  active  buffet  control 


Figure  26.  Effects  of  a closed- loop  active  buffet  control 
on  the  shock  position  - Experimental  results  (M=0.736; 
a = 2.5°;  dmean  = 15°;  Control  law:  A = 5;  cj)  = 90°; 
Cp(  § = 0.525)) 


Figure  28.  Effects  of  a closed- loop  active  buffet  control 
on  the  shock  position  - Direct  numerical  VIS  15  simula- 
tion (M=0.73;  a = 2.34°;  dmean  = 15°;  Control  law: 
A = 3;  0 = 45°;  Cp( f = 0.63)) 


Figure  27.  Effects  of  a closed-loop  active  buffet  control 
on  pressure  fluctuations  on  the  airfoil  - Experimental  re- 
sults (M=0.736;  a = 2.5°;  dmean  = 15°;  Control  law: 
A = 5;  0 = 30°;  Cp(f  = 0.525)) 
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Paper  #15 

Q by  Dr.  Bram  Elsenaar:  Is  the  preferred  frequency  of  75Hz  related  to  a preferred  structural  mode  of  the 
model,  like  bending  ? If  yes,  is  this  represented  in  the  theoretical  model;  or  can  you  explain  the  difference  in 
frequency  (50Hz  in  theory) 

A.  (Despre):  No  the  model  can  be  considered  as  stiff  (2-D  model  spanning  the  walls).  We  don’t  have  an 
explanation  for  this  difference;  it  may  be  related  to  the  difference  in  shock  position,  or  maybe  it  is  due  to  the 
TSP  inviscid  truncation  still  used  in  the  case  of  VIS  15,  or  to  the  form  of  TSP  variant  used. 
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